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The confusion and ambiguity encountered by students, in understanding virtual displacement 
and virtual work, is discussed in this article. A definition of virtual displacement is presented that 
allows one to express them explicitly for holonomic (velocity independent), non-holonomic (veloc- 
ity dependent), scleronomous (time independent) and rheonomous (time dependent) constraints. 
It is observed that for holonomic, scleronomous constraints, the virtual displacements are the dis- 
placements allowed by the constraints. However, this is not so for a general class of constraints. 
For simple physical systems, it is shown that, the work done by the constraint forces on virtual 
displacements is zero. This motivates Lagrange's extension of d'Alembert's principle to system of 
particles in constrained motion. However a similar zero work principle does not hold for the allowed 
displacements. It is also demonstrated that d'Alembert's principle of zero virtual work is necessary 
for the solvability of a constrained mechanical problem. We identify this special class of constraints, 
physically realized and solvable, as the ideal constraints. The concept of virtual displacement and 
the principle of zero virtual work by constraint forces are central to both Lagrange's method of 
undetermined multipliers, and Lagrange's equations in generalized coordinates. 
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I. INTRODUCTION 

Almost all graduate level courses in classi- 
cal mechanics include a discussion of virtual 
displacementiiSiiiiSiSiLSiSiiSiii and Lagrangian 
dynamico^i^i'^i'^i'^i^i'^i^i^i^°i"i^^. From the concept of 
zero work by constraint forces on virtual displacement, 
the Lagrange's equations of motion are derived. 

However, the definition presented in most accessible 
texts often seem vague and ambiguous to students. Even 
after studying the so called definition, it is rather com- 
monplace that a student fails to identify, whether a sup- 
plied vector is suitable as a virtual displacement, for a 
given constrained system. Though some of the more ad- 
vanced and rigorous treatisoi^ii present a more precise 
and satisfactory treatment, they are often not easily com- 
prehensible to most students. In this article we attempt 
a simple, systematic and precise definition of virtual dis- 
placement, which clearly shows the connection between 
the constraints and the corresponding allowed and virtual 
displacements. This definition allows one to understand 
how far the virtual displacement is 'arbitrary' and how 
far it is 'restricted' by the constraint condition. 

There are two common logical pathways of arriving at 
Lagrange's equation. 

1. Bernoulli's principle of virtual velocitjii (1717), 
d'Alembert's principle of zero virtual workii^ 
(1743), Lagrange's generalization of d'Alembert's 
principle to constrained system of moving particles, 
and Lagrange's equations of motion (1788)iiS4ii. 



2. Hamilton's principle of least action£*i£ (1834), and 
variational approach to Lagrange's equation. 

The two methods are logically and mathematically in- 
dependent and individually self contained. The first 
method was historically proposed half a century earlier, 
and it presents the motivation of introducing the La- 
grangian as a new physical quantity. The second method 
starts with the Lagrangian and the related action as 
quantities axiomatically describing the dynamics of the 
system. This method is applied without ambiguity in 
some textaiSii^ and courses^^. However one also finds 
intermixing of the two approaches in the literature and 
popular texts, often leading to circular definition and re- 
lated confusion. A rational treatment demands an in- 
dependent presentation of the two methods, and then a 
demonstration of their interconnection. In the present 
article we confine ourselves to the first method. 

In this approach, due to Bernoulli, d'Alembert and La- 
grange, one begins with a constrained system, defined by 
equations of constraints connecting positions, time and 
often velocities of the particles under consideration. The 
concept of virtual displacement is introduced in terms 
of the constraint equations. The external forces alone 
cannot maintain the constrained motion. This requires 
the introduction of forces of constraints. The imposition 
of the principle of zero virtual work by constraint forces 
gives us a 'special class of systems', that are solvable. 

A proper definition of virtual displacement is necessary 
to make the said approach logically satisfactory. However 
the various definitions found in popular texts are often 
incomplete and contradictory with one another. These 
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ambiguities will be discussed in detail in the next section. 

In the literature, e.g., Greenwood^ Eq.1.26 and Parsi^ 
Eq. 1.6.1, one encounters holonomic constraints of the 
form: 

(j>.j{xi,X2,.-.,X3N,t) ^0, j = l,2,...,k. (1) 

The dilferential form of the above equations are sat- 
isfied by allowed infinitesimal displacements {dxi}, 
(Greenwood^ Eq.1.27; Parsii Eq.1.6.3). 

^1^'^' + = ^ = 1,2,...,^. (2) 

For a system under above constraints the virtual displace- 
ments {Sxi}, satisfy the following equations (Greenwood^ 
Eq.1.28; Parsii Eq.1.6.5), 

3N „ , 

E|^'5^* = 0' J = 1,2,..., A;. (3) 

2—1 

The differential equations satisfied by allowed and virtual 
displacements arc different even for the non-holonomic 
case. Here, the equations satisfied by the allowed dis- 
placements {dxi} are (Goldsteini Eq.2.20, Greenwood^ 
Eq.1.29 and Parsii Eq.1.7.1), 

3JV 

ajidxi + ajtdt = , j = 1, 2, . . . , m. (4) 

i=l 

Whereas, the virtual displacements {Sxi} satisfy 
(Goldsteini Eq.2.21, Greenwood^ Eq.1.30; Parsi^ 
Eq.1.7.2), 

3N 

ajiSxi = , j = 1, 2, . . . , m. (5) 

i=l 

Thus, there appear in the literature certain equations, 
namely Eq.© and Eq.Q, which are always satisfied by 
the not so precisely defined virtual displacements. It may 
be noted that these equations are connected to the con- 
straints but are not simply the infinitesimal forms of the 
constraint equations, i.e., Eqs.Q and This fact is 
well documented in the literatureitSii^. However, the na- 
ture of the difference between these sets of equations, 
i.e, Eqs.(|3J) and © on one hand and Eqs.Q and on 
the other, and their underlying connection are not ex- 
plained in most discussions. One may consider Eq.Q 
or Eq.Q, as independent defining equation for virtual 
displacement. But it remains unclear as to how, the vir- 
tual displacements {Sxi} defined by two different sets of 
equations for the holonomic and the non-holonomic cases, 
viz., Eq.|(21) and Eq.jSJ), correspond to the same concept 
of virtual displacement. 

We try to give a physical connection between the defi- 
nitions of allowed and virtual displacements for any given 
set of constraints. The proposed definition of virtual dis- 
placement (Sec. HA) as difference of two unequal allowed 



displacements (satisfying Eq.® or Eq.Q) over the same 
time interval; automatically ensures that virtual displace- 
ments satisfy Eq.© and Eq.® for holonomic and non- 
holonomic systems respectively. We show that in a num- 
ber of natural systems, e.g., pendulum with fixed or mov- 
ing support, particle sliding along stationary or moving 
frictionless inclined plane, the work done by the forces 
of constraint on virtual displacements is zero. We also 
demonstrate that this condition is necessary for the solv- 
ability of a constrained mechanical problem. Such sys- 
tems form an important class of natural systems. 

A. Ambiguity in virtual displacement 

In the literature certain statements appear in reference 
to virtual displacement, which seem confusing and mu- 
tually inconsistent, particularly to a student. In the fol- 
lowing we present few such statements found in common 
texts. 

1. It is claimed that (i)a virtual displacement Sr is 
consistent with the forces and constraints imposed 
on the system at a given instant ii; (ii) a virtual dis- 
placement is an arbitrary, instantaneous, infinites- 
imal change of position of the system compatible 
with the conditions of constraint; (iii) virtual dis- 
placements are, by definition, arbitrary displace- 
ments of the components of the system, satisfy- 
ing the constraint; (iv) virtual displacement does 
not violate the constraintsi^; (v) we define a vir- 
tual displacement as one which does not violate the 
kinematic relations^; (vi) the virtual displacements 
obey the constraint on the motion^. These state- 
ments imply that the virtual displacements satisfy 
the constraint conditions, i.e., the constraint equa- 
tions. However this is true only for holonomic, scle- 
renomous constraints. We shall show that for non- 
holonomic constraints, or rheonomous constraints, 
e.g., a pendulum with moving support, this defini- 
tion violates the zero virtual work principle. 

2. It is also stated that (i) virtual displacements do 
not necessarily conform to the constraints^; (ii) the 
virtual displacements Sq have nothing to do with 
actual motion. They are introduced, so to speak, as 
test quantities, whose function it is to make the sys- 
tem reveal something about its internal connections 
and about the forces acting on it; (iii) the word 
"virtual" is used to signify that the displacements 
are arbitrary, in the sense that they need not corre- 
spond to any actual motion executed by the systent; 
(iv) it is not necessary that it (virtual displacement) 
represents any actual motion of the system^; (v) 
it is not intended to say that such a displacement 
(virtual) occurs during the motion of the particle 
considered, or even that it could occuf^; (vi) virtual 
displacement is any arbitrary infinitesimal displace- 
ment not necessarily along the constrained patlt. 
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From the above we understand that the virtual dis- 
placements do not neeessarily satisfy the constraint 
equations, and they need not be the ones actually 
realized. We shall see that these statements are 
consistent with physical situations, but they cannot 
serve as a satisfactory definition of virtual displace- 
ment. Statements like: "not necessarily conform to 
the constraints" or "not necessarily along the con- 
strained path" only tell us what virtual displace- 
ment is not, they do not tell us what it really is. 
Reader should note that there is a conflict between 
the statements quoted under items 1 and 2. 

Thus it is not clear from the above, whether the 
virtual displacements satisfy the constraints, i.e., 
the constraint equations, or they do not. 

3. It is also stated that (i) virtual displacement is to 
be distinguished from an actual displacement of the 
system occurring in a time interval dii; (ii) it is 
an arbitrary, instantaneous, change of position of 
the systeni; (iii) virtual displacement 5r takes place 
without any passage of timsi^. (iv) virtual displace- 
ment has no connection with the time - in contrast 
to a displacement which occurs during actual mo- 
tion, and which represents a portion of the actual 
patfA; (v) one of the requirements on acceptable vir- 
tual displacement is that the time is held fixe(&. We 
even notice equation like : "Jx^ = dxi for dt = 0"—. 
The above statements are puzzling to a student. If 
position is a continuous function of time, a change 
in position during zero time has to be zero. In 
other words, this definition implies that the vir- 
tual displacement cannot possibly be an infinitesi- 
mal (or differential) of any continuous function of 
time. In words of Arthur Haas: since its (virtual 
displacement) components are thus not functions of 
the time, we are not able to regard them as differ- 
entials, as we do for the components of the element 
of the actual pathm. It will be shown later (Sec. II), 
that virtual displacement can be looked upon as a 
differential. It is indeed a differential change in po- 
sition or an infinitesimal displacement, consistent 
with virtual velocity vj, (t) , taken over a time inter- 
val dt (see Eq.lHSJ). 

4. Virtual displacement is variously described as: ar- 
bitrary, virtual, and imaginary^i^i^i^. These ad- 
jectives make the definition more mysterious to a 
student. 

Together with the above ambiguities, students are of- 
ten unsure whether it is sufficient to discuss virtual dis- 
placement as an abstract concept or it is important to 
have a quantitative definition. Some students appreciate 
that the virtual displacement as a vector should not be 
ambiguous. The principle of zero virtual work is required 
to derive Lagrange's equations. For a particle under con- 
straint this means that the virtual displacement is always 
orthogonal to the force of constraint. 



At this stage a student gets further puzzled. Should he 
take the forces of constraint as supplied, and the principle 
of zero virtual work as a definition of virtual displacement 
? In that case the principle reduces merely to a defini- 
tion of a new concept, namely virtual displacement. Or 
should the virtual displacement be defined from the con- 
straint conditions independently ? The principle of zero 
virtual work may then be used to obtain the forces of 
constraint. These forces of constraint ensure that the 
constraint condition is maintained throughout the mo- 
tion. Hence it is natural to expect that they should be 
connected to and perhaps derivable from the constraint 
conditions. 



II. VIRTUAL DISPLACEMENT AND FORCES 
OF CONSTRAINT 

A. Constraints and Virtual displacement 

Let us consider a system of constraints that are ex- 
pressible as equations involving positions and time. They 
represent some geometric restrictions (holonomic) either 
independent of time (sclerenomous) or explicitly depen- 
dent on it (rheonomous). Hence for a system of N par- 
ticles moving in three dimensions, a system of (s) holo- 
nomic, rheonomous constraints are represented by func- 
tions of {rfc} and {t), 

/i(ri,r2, . . . ,rAr,t) = 0, i=l,2,...,s. (6) 

The system may also be subjected to non-holonomic con- 
straints which are represented by equations connecting 
velocities {v^}, positions {rfc} and time {t). 

N 

^ A,fc ■ Vfc + , i = 1, 2, . . . , m, (7) 

fc=i 

where {Aik} and {Ait} are functions of positions 
{ri, r2, . . . , rjv} and time (t). The equations for non- 
holonomic constraints impose restrictions on possible or 
allowed velocity vectors {vi, V2, . . . , vjv}, for given posi- 
tions {ri, r2, . . . , rjv} and time (t). The holonomic con- 
straints given by Eq.®, are equivalent to the following 
equations imposing further restriction on the possible or 
allowed velocities. 

For a system of N particles under (s) holonomic 
and (to) non-holonomic constraints, a set of vectors 
{vi, V2, . . . , Vat} satisfying Eq.Q and Eq.Q are called 
allowed velocities. It is worth noting at this stage that 
there are many, in fact infinitely many, allowed velocities, 
since we have imposed only {s-\-m) number of scalar con- 
straints, Eq.([7I) and Eq.®, on (3A^) scalar components 
of the allowed velocity vectors. 
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At any given instant of time, the difference of any 
two such non-identical allowed sets of velocities, inde- 
pendently satisfying the constraint conditions, are called 
virtual velocities. 



1,2,. 



An infinitesimal displacement over time [t, t + dt), due to 
allowed velocities, will be called the allowed infinitesimal 
displacement or simply allowed displacement. 



drk = ^k{t)dt, 



fc = l,2. 



,N. (9) 



Allowed displacements {dr^} together with differential of 
time [dt) satisfy the infinitesimal form of the constraint 
equations. From Eq.® and Eq.©, we obtain. 



^ / o J- 



fc=l ^ ^ 

N 

^ Ajfc • dvk + A^tdt = , 
fe=i 



z= 1,2,. 



,s, (10) 



i 1,2, . . . ,m. (11) 



holonomic 









( = 1,2,...,.? 



constraints 



non-holonomic 

N 

y A. -v +A =0 

^-^ ik k It 

k=l 

i = l ,2 , m 



df df 

I 1 , I 



t=i 9r, 



dr, + — -dt = 
dt 

(=1,2,...,* 



Eqns X ' 

for X^,-dr^+A_^dt=0 
allowed k=i 
displ. 



i=l ,2 , ... ,m 



5r, =dr, —dr' 

k k k 



y '_ |5 J. _ Q ^ Equations 

' virtual 
(=1,2, .... * displacements 



k = \ 

i = \ ,2 , m 



FIG. 1: Virtual displacement defined as difference of allowed 
displacements 



As there are many independent sets of allowed veloci- 
ties, we have many allowed sets of infinitesimal displace- 
ments. We propose to define virtual displacement as the 
difference between any two such (unequal) allowed dis- 
placements taken over the same time interval (t,t + dt), 



6vk = drk - drL 



= 1,2 



, ij, . . . , 



TV. (12) 



Thus virtual displacements are infinitesimal displace- 
ments over time interval dt due to virtual velocity Vfe(i), 



The absence of the {dfi / dt) and An in the above equa- 
tions, Ea. H14() and Ea. (|15() . gives the precise meaning to 
the statement: ^^virtual displacements are the allowed dis- 
placements in the case of frozen constraints" . The con- 
straints are frozen in time in the sense that we make 
the (dfi/dt) and terms zero, though the dfi/dri. 
and Aik terms still involve both position {ri, r2, . . . , fat}, 
and time (t). In the case of stationary constraints, i.e., 
/i(ri, . . . ,rAr) = 0, and J2k-^ik{'^ii ■ ■ • ii"Ar) • Vfc = 0, the 
virtual displacements are identical with allowed displace- 
ments as {dfi/dt) and Ait are identically zero. 



5rk^vk{t)dt^ {vk{t)-Vkit))dt, fc = l,2,...,iV. (13) 

This definition is motivated by the possibility of (i) iden- 
tifying a special class of ^ ideal constraints'" (Sec.IIC), and 
(ii) verifying ' the principle of zero virtual worK in com- 
mon physical examples (Sec. III). It may be noted that, 
by this definition, virtual displacements {5vk} are not in- 
stantaneous changes in position in zero time. They arc 
rather smooth, differentiable objects. 

The virtual displacements thus defined, satisfy the ho- 
mogeneous part of the constraint equations, i.e., Ea. l(Tn|) 
and Eq.imj with dfi/dt = and Au = 0. Hence, 



N 



> — • 5Tk = 0, 

N 

^ A,k ■ 5vk = 0, 



z = l,2,...,s, (14) 



i = l,2,...,m. (15) 



fe=i 



The logical connection between the equations of con- 
straint, equations for allowed displacements and equa- 
tions for virtual displacements are presented in FIG^ 



B. Existence of forces of constraints 

In the case of an unconstrained system of N particles 
described by position vectors {rfc} and velocity vectors 
{vfe}, the motion is governed by Newton's Law, 



TOfcafc ^ Ffc(r;,V(,t), 



fc,Z= 1,2,...,7V (16) 



where mk is the mass of the k*'* particle, is its accel- 
eration and Ffc is the total external force acting on it. 
However, for a constrained system, the equations of con- 
straint, namely Eq.® and Eq.lQ, impose the following 
restrictions on the allowed accelerations. 



^a/, ^dfdfA ^d(df,. ^ 



z= 1,2,, 



N 



Ed . d . ^ 

A^fc ■ ELk + — A,,fc ■ Vfc + —Ait = 0, 

k=l 



dt 



dt 



* = 1,2,, 



,s (17) 



(18) 
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Given {r^}, {v^} one is no longer free to ehoose all the 
accelerations {sik} independently. Therefore in general 
the accelerations {aj.} allowed by Ea. H17|) . Ea. ll8|) are 
incompatible with Newton's Law, i.e., Ea. (|16|) . 

This implies that during the motion the constraint con- 
dition cannot be maintained by the external forces alone. 
Physically some additional forces, e.g., normal reaction 
from the surface of constraint, tension in the pendulum 
string, come into play to ensure that the constraints are 
satisfied throughout the motion. Hence one is compelled 
to introduce forces of constraints {R/c} and modify the 
equations of motion as. 



1,2,..., TV. (19) 



force is zero, on allowed as well as virtual displacement. 



N 



N 



^Rfc-drfc = 0, ^Rfc-<5r, 



0. 



k=l 



k=l 



When the constraint surface is in motion, the allowed 
velocities, and hence the allowed displacements are no 
longer tangent to the surface (see Sec. III). The virtual 
displacement however remains tangent to the constraint 
surface. As the surface is frictionless, it is natural to 
assume that the force of constraint is still normal to the 
instantaneous position of the surface. Hence the work 
done by normal reaction on virtual displacement is zero. 
However the work done by constraint force on allowed 
displacements is no longer zero. 



Newton's equation of motion 
m,r=F,, k = l,2,- N 

k k k 



forces of constraint 
remove the 
incompatibility 



Equations of 
constraint 



N 



N 



Oj 









i=\,2,...,s 



Sa -v +A. =0 

(= 1 , 2, ... , m 



Newton's equation of motion 
with constraint forces. 
'"*'^, = F^ + R,, k=xx..N 



FIG. 2: Existence of constraint forces 

FIG [3 presents the connection between the equations 
of constraints and forces of constraints. 

Now the problem is to determine the motion of N par- 
ticles, namely their positions {rk{t)}, velocities {vfc(t)} 
and the forces of constraints {Rfc}, for a given set of 
external forces {F^}, constraint equations, Eq.JHll and 
Eq.©, and initial conditions {rfc(O), v;;(0)}. It is im- 
portant that the initial conditions are also compatible 
with the constraints. There are a total of (6iV) scalar 
unknowns, namely the components of Yk{t) and R^, con- 
nected by (3A^) scalar equations of motion, Ea. (|19|) . and 
(s + m) equations of constraints, Eq.® and Eq.Q. For 
(6iV > 3A^ -l- s -|- m) we have an under-determined sys- 
tem. Hence to solve this problem we need (3N — s — m) 
additional scalar relations. 



C. Solvability and ideal constraints 

In simple problems with stationary constraints, e.g., 
motion of a particle on a smooth stationary surface, we 
observe that the allowed displacements are tangential to 
the surface. The virtual displacement being a difference 
of two such allowed displacements, is also a vector tan- 
gential to it. For a frictionless surface, the force of con- 
straint, the so called 'normal reaction', is perpendicular 
to the surface. Hence the work done by the constraint 



^Rfe-drfc^O, ^Rfc-(5rfe = 0. 



(20) 



k=l 



fc=l 



In a number of physically interesting simple problems, 
such as, motion of a pendulum with fixed or moving sup- 
port, motion of a particle along a stationary and moving 
slope, we observe that the above interesting relation be- 
tween the force of constraint and virtual displacement 
holds (see Sec. HI). As the 3iV scalar components of the 
virtual displacements {(Jr^} are connected by (s + m) 
equations, Ea. ((n|l and Ea. lfT^ . only n = {3N — s — m) of 
these scalar components are independent. If the (s -I- m) 
dependent quantities are expressed in terms of remaining 
(3A^ — s — m) independent objects we get. 



0. 



(21) 



where {xj} are the independent components of {rfc}. 
{Rj} are the coefficients of {Sxj}, and are composed of 
different {Rfc}. Since the above components of virtual 
displacements {6xj} are independent, one can equate 

each of their coefficients to zero {Rj = 0). This brings in 
exactly {3N — s — m) new scalar conditions or equations 
that are needed to make the system solvable (see FIG|2l. 

Thus we have found a special class of constraints, 
which is observed in nature (Sec. HI) and which gives us a 
solvable mechanical system. We call this special class of 
constraints, where the forces of constraint do zero work 
on virtual displacement, i.e., Rfc • Jrfc = 0, the ideal 
constraint. 

Our interpretation of the principle of zero virtual work, 
as a definition of an ideal class of constraints, agrees with 
Sommerfcld. In his exact words, "a genera? postulate of 
mechanics: in any mechanical systems the virtual work 
of the reactions equals zero. Far he it from us to want to 
give a general proof of this postulate, rather we regard it 
practically as definition of a mechanical system"—. 
(Boldface is added by the authors). 
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Satisfied by simple ptiysical 
systems (Sec III) 



Equations for virtual 
displacements 5r 

^ k 

k=-i or, 



Newton's equation of motion 
with constraint forces 
m,f =F,+R,, k = l,...N 

k k k 





Together 


( — 


form a 




solvable 




problem 



Satisfied by simple physical 
systems (Sec III) 

d' Alembert's principle 
(Lagrange's extension) 

ER^.5r^=0 

k=l 

for system of moving particles 



FIG. 3: Solvability under ideal constraints 

III. EXAMPLES OF VIRTUAL 
DISPLACEMENTS 

A. Simple Pendulum with stationary support 

The motion of a pendulum is confined to a plane and 
its bob moves keeping a fixed distance from the point 
of suspension (see FIGQ. The equation of constraint 
therefore is, 

/(x, t) - a;2 + y2 _ ^2 ^ 0, 
where tq is the length of the pendulum. Whence 
df 9J df 



dx 



2x, 



dy 



22/, 



dt 



0. 



The constraint equation for allowed velocities, Eq.®, 
becomes, 

X -Vx+y -Vy =0. 

Hence the allowed velocity {vx, Vy) is orthogonal to the 
instantaneous position (x^y) of the bob relative to sta- 
tionary support. The same may also be verified taking a 
plane polar coordinate. 

The allowed velocities and the allowed displacements 
are perpendicular to the line of suspension. The virtual 
velocities and the virtual displacements, being the differ- 
ence of two unequal allowed velocities and displacements 
respectively, are also perpendicular to the line of suspen- 
sion. 

dr ^ v{t)dt, dr' = v'{t)dt, 

Sr = (v(t) - v'{t))dt. 




/P' dr'=v'(t)dt 
^ dr=v(t)dt 



FIG. 4: Allowed and virtual displacements of a pendulum 
with stationary support 



Although the virtual displacement is not uniquely speci- 
fied by the constraint, it is restricted to be in a plane per- 
pendictilar to the instantaneous line of suspension. Hence 
it is not 'completely arbitrary'. 

The ideal string of the pendulum provides a tension 
(T) along its length, but no shear. The work done by 
this tension on botlr allowed and virtual displacements is 
zero, 



T • dr = 0, 



T • (5r = 0. 



B. Simple Pendulum with moving support 

Let us first consider the case when the support is mov- 
ing vertically with a velocity u. The motion of the pen- 
dulum is still confined to a plane. The bob moves keep- 
ing a fixed distance from the moving point of suspension 
(FIGEI). The equation of constraint is, 

!{x,y,t) ^ x" + {y - utf - rl=^, 

where it is the velocity of the point of suspension along 
a vertical direction. 
Whence 

1 = 2., |=2(,-.t), %^-2ui^y-ut). 

Hence the constraint equation gives, 
X ■ Vx + {y — ut) ■ Vy — u{y — ut) = 0, 
or, 

X ■ Vx + {y — ut) ■ (Vy — u) = 0. 

The allowed velocities {vx, Vy) and the allowed displace- 
ments, are not orthogonal to the instantaneous position 
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FIG. 5: Allowed and virtual displacements of a pendulum 
with moving support 

of the bob relative to the instantaneous point of suspen- 
sion (x, y — ut). It is easy to verify from the above equa- 
tion that the allowed velocity [vx, Vy) is equal to the 
sum of a velocity vector (vx, Vy — u) perpendicular to 
the position of the bob relative to the point of suspen- 
sion {x,y — ut), and the velocity of the support (0, u). 
If we denote v{t) = {vx,Vy), v±{t) = {vx,Vy — u) and 
u — (0, u), then, 

V(t) = Vi(i)+U. 

The allowed displacements are vectors collinear to al- 
lowed velocities. A virtual displacement being the differ- 
ence of two allowed displacements, is a vector collinear to 
the difference of allowed velocities. Hence it is orthogonal 
to the instantaneous line of suspension. 

dr = v{t)dt = w^{t)dt + udt, 
dr' = V{t)dt = v'^{t)dt + udt, 
Sr = (v(t) - v'{t))dt = (vi(<) - v'^{t))dt. 

Hence none of these allowed or virtual vectors are 'arbi- 
trary'. 

At any given instant, an ideal string provides a ten- 
sion along its length, with no shear. Hence the constraint 
force, namely tension T, does zero work on virtual dis- 
placement. 

T ■ dr = T -udt ^0, T ■ Sr = T ■ v± dt ^ 

For the support moving in a horizontal or any arbitrary 
direction, one can show that the allowed displacement 
is not normal to the instantaneous line of suspension. 
But the virtual displacement, as defined in this article, 
always remains perpendicular to the instantaneous line 
of support. 



C. Motion along a stationary inclined plane 

Let us consider a particle sliding along a stationary 
inclined plane as shown in FIG El The constraint here is 
more conveniently expressed in polar coordinates. The 
constraint equation is, 

fir,e) = 6-00 = 

where Oq is the angle of the slope. Hence the constraint 
equation for allowed velocities, Eq.®, gives, 

fdf\ ^df fdf\ fdf\ df 

= • r + 1 • {r6) +0 = 

Hence 9 = Q, implying that the allowed velocities are 
along the constant 9 plane. Allowed velocity, allowed 
and virtual displacements are, 

V = rf, dr = frdt, 

v' = r'r, dr' = f'rdt, 

6r = {v — v')dt = (r — r')r dt. 

where f is a unit vector along the slope. 




FIG. 6: Allowed and virtual displacements of a particle sliding 
along a stationary slope 

As the inclined slope is frictionless, the constraint force 
N is normal to the surface. The work done by this force 
on allowed as well as virtual displacement is zero. 

N • dr = 0, N • (5r = 

D. Motion along a moving inclined plane 

For an inclined plane moving along the horizontal side 
(FIG|7I), the constraint is given by, 

[x — ut) , 

- cot(6lo) = 0, 

V 

f{x,y) = {x ~ ut) - cot{9o)y = 0. 

Whence the constraint for allowed velocities Eq.® be- 
comes, 

{x — u) — cot{9o)y = 0. 
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<5r 



d r^v'dt 



V— V 



v'„dt\ 



= vdt 



v„dr 



FIG. 7: Allowed and virtual displacements of a particle sliding 
along a moving slope 



Hence the allowed velocity (i, y) is the sum of two 
vectors, one along the plane {i — u,y), and the other 
equal to the velocity of the plane itself {u, 0). 

V(i) = V||(t)+U. 

Allowed displacements are vectors along the allowed ve- 
locities, however the virtual displacement is still a vector 
along the instantaneous position of the plane. 

dr = (v|| (t) + u)dt, dr' = (vj, [t) + u)dt, 

Sr = (v(t)-v'(0)dt = (v||(0-v||(i))di. 

For the moving frictionless slope, the constraint force pro- 
vided by the surface is perpendicular to the plane. Hence 
the work done by the constraint force on virtual displace- 
ment remains zero. 



N • dr ^ 0, 



N • = 0. 



IV. LAGRANGE'S METHOD OF 
UNDETERMINED MULTIPLIERS 

A constrained system of particles obey the equations 
of motion given by. 



rukSik = Ffc + Rfc 



1,2,. 



where is the mass of the k**^ particle, is its accel- 
eration. Ffc and Rfc arc the total external force and force 
of constraint on the k*'' particle. If the constraints are 
ideal, we can write. 



N 

fc=i 

whence wc obtain, 

N 

^(mfcafc - Ffe) • Srk = 0. 

fc=i 



(22) 



(23) 



If the components of {Svk} were independent, we could 
recover Newton's Law for unconstrained system from this 



equation. However for a constrained system {Svk} are de- 
pendent through the constraint equations, Ea. (|14() and 
Ea. (|15|) . for holonomic and non-holonomic systems re- 
spectively. 



N 



k=l 
N 



z = l,2,. 



(14) 



(15) 



k=l 



We multiply Ea. lfTH) successively by (s) scalar multi- 
pliers {Ai, A2, . . . As}, Ea. ((T5|l successively by (m) scalar 
multipliers {/ii, ^2, ■ ■ ■ Mm} ^-^d then subtract them from 
the zero virtual work equation, namely Ea. (|22|) . 



N 



k=l 



(24) 



These multipliers {A^} and {/ij} are called the Lagrange's 
multipliers. Explicitly in terms of components. 



N 

E 



i=l 



dxk 



i=i 



6xk 



N 



N 



-Y^lYhSyk + Y.[Z]k5zk = 0, 



(25) 



fc=i 



fe=i 



where [Y]k and [Z]k denote the coefficients of Syk and 
Szk respectively. 

The constraint equations, Ea. (|14|l and Eg. ljlSII . allow 
us to write the (s -f- m) dependent virtual displacements 
in terms of the remaining n = {3N — s — to) independent 
ones. We choose (s -|- to) multiphers {Ai, A2, . . . , A^} and 
{/^i: A*2, • ■ • , Mm}, such that the coefficients of (s + to) 
dependent components of virtual displacement vanish. 
The remaining virtual displacements being independent, 
their coefficients must vanish as well. Thus it is possible 
to choose {Ai, A2, . . . , A^} and {^i, ^J■2, ■ • ■ , ^J■m} such that 
all coefficients {[X]fc,[y]fc,[Z]fc} of virtual displacements 
{dxk,Syk,Szk} in Ea. (|25|) vanish. Hence we can express 
the forces of constraint in terms of the Lagrange's multi- 
pliers. 



Rfc = ^ A, 



dvk 



1,2,, 



,N (26) 



Thus the problem reduces to finding a solution for the 
equations. 



d f- 

iTikSik = Ffc+^ A,^+^yUiAjfc,A: = 1,2,, 



together with the equations of constraint, 
/,(ri,r2, . . . ,rAr,t) = 0, 



1,2,, 



,N (27) 



,,5 (6) 
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Eq.s for constrained system 

N 



lc=\ 
N 



N * 



yA,,5r=0, i = \.2,....m 



Lagrange's equations 
(method of undetermined multipliers) 
■AS/™ 

,= 1 or 



k = \,l,...,N 



1=1 



/,(r,,r,,...,r^,0 = 0, 

N 

y A -v +A. =0, 

ik k It ' 



i = \.2.....s 
i = \ ,2 , .. ,m 



No. of unknowns 
= No. of Eq.s 



Unknown 


No. of scalar 


quantities 


components 






A 






m 



FIG. 8: Lagrange's method of undetermined multipliers; solv- 
ability 



and 

N 

^ A,fc • vfc + A,t = , 
fc=i 



z = l,2,. 



(7) 



Here we have to solve (3iV + s + m) scalar equations 
involving (iN +s+m) unknown scalar quantities, namely 
{xk,yk, Zk, K, fJ'j} (see FIG|Sll. After solving this system 
of equations for {xk,yki Zk, Ki ^J■j}, one can obtain the 
forces of constraint {Rfc} using Ea. l(^ . 



V. LAGRANGE'S EQUATIONS IN 
GENERALIZED COORDINATES 

For the sake of completeness we discuss very briefly 
Lagrange's equations in generalized coordinates. A 
more complete discussion can be found in most 
textoii^iVi5iW,9ii°i"ii2ii3ii4^ Consider a system of N 

particles under (s) holonomic and (to) non-holonomic 
constraints. In certain suitable cases, one can express 
(s + to) dependent coordinates in terms of the remaining 
(3iV — s — to) independent ones. It may be noted that 
such a complete reduction is not possible for general cases 
of non-holonomic and time dependent constraintaiil. If 
we restrict our discussion to cases where this reduction is 
possible, one may express all the 3N scalar components 



of position {ri, r2, . . . , r^r} in terms of (3A^ — s — to) in- 
dependent parameters {gi, (72, • • • , 9n} and time (t). 

rk ^ rk{qi,q2, ■ ■ ■ ,qn,t), k = 1,2, . . . ,N (28) 

The allowed and virtual displacements are given by, 

" 9rfc , , dvk 



dv. 



dt 



^ % " 

srk = y^^^ij. 



From Ea. (|23|l we obtain, 
dtk 



-dt. 



k = 1,2,. ..,N 



k^l,2,...,N. (29) 



N 



fc=l 



N 



dvk 



k=l 



(30) 



Introduce the expression of kinetic energy. 



1 ^ 

T=\Y^^kil, 

fc=i 



and that of the generalized force, 

N 



dvk 



j = l,2,...,n. (31) 



k=i 

After some simple algebra one finds. 



^(ddT dT . 



(32) 



As {qi,q2, ■ ■ ■ ,qn} are independent coordinates, coeffi- 
cient of each 5qj must be zero separately. Hence (see 
FIG®, 



d dT dT _ ^ 
dt dqj dqj ' 



j = l,2,...,n. (33) 



These are called Lagrange's equations in generalized co- 
ordinates. To proceed further one has to impose addi- 
tional conditions on the nature of forces {Fk} or {Qj}. 



In problems where forces {Fk} are derivable from a 
scalar potential y(ri,r2, . . . ,rAr), 

Ffe = -VfcF(ri,r2,...,rjv), k = l,2,...,N. (34) 

One can obtain the generalized force as, 
'drk\ dV 



Qj = -VkV 



dqj 



dqj' 



j = l,2,...,n. (35) 



10 



Constrained system (holonomic) 

N 

Z(«.f.-Fj.5r =0 



«:= 1 

" df 



X— 5r^ = 0, , = 1,2,.... 



virtual displacements in terms of 
generalized coordinates. 

6r = > 6q , i = 'i,...,N 

i da 



Lagrange's equation in generalized coordinates 


d 


dT 




8T 












dt 






dq 



FIG. 9: Lagrange's equations in generalized coordinates 



Where V is the potential V expressed as a function of 
92, • ■ • , ^n}- In addition as the potential V is in- 
dependent of the generalized velocities, we obtain from 
Eq.lEHJl, 



dd{T-V) d{T-V) 



dt dqj 



0, j = l,2,...,n. (36) 



At this stage one introduces the Lagrangian function, 
L — T ~V . In terms of the Lagrangian, the equations of 
motion take up the form. 



d dL dL _^ 
dt dqj dqj ' 



j = l,2,...,n. (37) 



VI. CONCLUSION 

Iir this article we make an attempt to present a quan- 
titative definition of the virtual displacement. We show 
that for certain simple cases the work done by the forces 
of constraint on virtual displacement is zero. We also 
demonstrate that this zero virtual work principle gives 
us a solvable class of problems. Hence we define this 
special class of constraint as the ideal constraint. We 
demonstrate in brief how one can solve a general mechan- 
ical problem by: i) Lagrange's method of undetermined 
multipliers and ii) Lagrange's equations in generalized 
coordinates. 

In the usual presentations of Lagrange's equation 
based on virtual displacement and d'Alembert's princi- 
ple, Eq.(3) and Eq.(5) are satisfied by the virtual dis- 
placements. One may consider these equations as the 
definition of virtual displacements. However the situ- 
ation is far from satisfactory, as separate defining equa- 
tions are required for different classes of constraints. This 
adhoc definition also fails to clarify the actual connection 
between the virtual displacements and the equations of 
constraints. 



At this stage one introduces d'Alembert's principle of 
zero virtual work. Bernoulli^ (1717) and d'AlembertLi^ 
(1743) originally proposed this principle for a system in 
static equilibrium. The principle states that the forces of 
constraint do zero work on virtual displacement. For sys- 
tems in static equilibrium, virtual displacement meant an 
imaginary displacement of the system that keeps its stat- 
ical equilibrium unchanged. Lagrange generalized this 
principle to a constrained system of particles in motion. 
This principle is crucial in arriving at Lagrange's equa- 
tion. However, most texts do not clearly address the 
questions, (i) why one needs to extend d'Alembert's prin- 
ciple to particles in motion, and (ii) why the work done 
by constraint forces on virtual displacements, and not on 
allowed displacements, is zero ? 

In the present article the allowed infinitesimal displace- 
ments are defined as ones that satisfy the infinitesimal 
form of the constraint equations. They arc the displace- 
ments that could have been possible if only the con- 
straints were present. Actual dynamics, under the given 
external forces, would choose one of these various sets of 
displacements as actual displacement of the system. The 
definition of virtual displacement as difference of two un- 
equal allowed displacements over the same infinitesimal 
time interval {t, t+dt), gives a unified definition of virtual 
displacement. This definition of virtual displacement sat- 
isfies the appropriate equations found in the literature, 
for both holonomic and non- holonomic systems. 

It is shown that Newton's equation of motion with ex- 
ternal forces alone, is inconsistent with equations of con- 
straint. Hence the forces of constraint are introduced. 
Now there are 3A'^ equations of motion and (s 4- m) 
equations of constraint involving 6A^ unknown scalars 
{rfc(t),Rfe}. Without additional condition (d'Alembert 
principle), the problem is underspecified and unsolvable. 

It is verified that for simple physical systems, the 
virtual displacements, as defined in this article, sat- 
isfy d'Alembert principle for particles in motion. The 
rheonomous examples discussed in Sec. Ill show why the 
forces of constraint do zero work on virtual displace- 
ments and not on allowed displacements. The additional 
equations introduced by d'Alembert principle make the 
problem solvable. These justify (i) the peculiar defini- 
tion of virtual displacements and the equations they sat- 
isfy, (ii) Lagrange's extension of d'Alembert's principle 
to particles in motion, (iii) why the zero work princi- 
ple is related to virtual displacement and not to allowed 
displacemeirt. Once the system is solvable, two meth- 
ods, originally proposed by Lagrange can be used, and 
are demonstrated. For Lagrange's method of undeter- 
mined multipliers, one solves (3iV -I- s) equations to ob- 
tain the motion of the system {rk{t), . . . ,ri\[{t)} and La- 
grange's multipliers {Xi, . . . , As}. The forces of constraint 
{Rfc, . . . ,Riv} are expressed in terms of these multipli- 
ers. For Lagrange's equations in generalized coordinates 
one solves (3iV — s) equations to obtain the time evo- 
lution of the generalized coordinates {qj — qj{t), j = 
1, . . . ,3A^ ^ s}. This gives the complete description of 
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the motion {rfc = rk{qi,q2, .. .,qn,t), k ^ 1,. . .,N}, ig- 
noring the calculation of the constraint forces. It may be 
noted that about a century later AppcU's cquation3£ii^iS£ 
were introduced for efficiently solving non-holonomic sys- 
tems. 

It is interesting to note that both the above mentioned 
methods require the principle of zero virtual work by con- 
straint forces as a crucial starting point. In the case of 
Lagrange's method of undetermined multipliers we start 
with the ideal constraint condition Ea. (|22|l . From there 
we obtain Ea. (P^ -Ea. lf77|) . Ea. (PS|) expresses the con- 
straint forces in terms of Lagrange's multipliers. For La- 
grange's equations in generalized coordinates, we start 
with the ideal constraint, Ea. H22|) . Wc work our way 
through Eq.^, Eq.ljSOJ), Eq.^ and finally obtain La- 
grange's equations in generalized coordinates, Eq.|j22I) 
and Eq.^^. The last figure, FIG.lUni), gives the com- 
plete logical flow of this article. 
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Newton's equation of motion 
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Newton's equation of motion 
with constraint forces. 
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Equations for 
virtual displacements 



= 0, 



;= 1,2,. .,s 



A/ Together \A 
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FIG. 10: Logical connection between constraint equations, virtual displacements, forces of constraints, d' Alembert's principle 
and Lagrange's equations. 



